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ABSTRACT 

We calculate by means of the Press-Schechter formalism the density profile developed by dark-matter 
halos during accretion, i.e., the continuous aggregation of small clumps. We find that the shape of the 
predicted profile is similar to that shown by halos in high-resolution cosmological simulations. Fur- 
thermore, the mass-concentration relation is correctly reproduced at any redshift in all the hierarchical 
cosmologies analyzed, except for very large halo masses. The role of major mergers, which can cause the 
rearrangement of the halo structure through violent relaxation, is also investigated. We show that, as a 
result of the boundary conditions imposed by the matter continuously infalling into the halo during the 
violent relaxation process, the shape of the density profile emerging from major mergers is essentially 
identical to the shape the halo would have developed through pure accretion. This result explains why, 
according to high-resolution cosmological simulations, relaxed halos of a given mass have the same den- 
sity profile regardless of whether they have had a recent merger or not, and why both spherical infall 
and hierarchical assembly lead to very similar density profiles. Finally, we demonstrate that the density 
profile of relaxed halos is not affected either by the capture of clumps of intermediate mass. 

Subject headings: cosmology: theory — dark matter — galaxies: halos, structure, formation 



1. INTRODUCTION 

The dominant dark component of matter in the universe 
is clustered in bound halos which form the skeleton of all 
astronomical objects of cosmological interest, from dwarf 
galaxies to rich galaxy clusters. The determination of the 
inner structure of such halos has been addressed both an- 
alytically and by numerical simulations. 

The analytical approach to this problem was pioneered 
by Gunn & Gott (1972), who adopted the simphfying as- 
sumption that halos grow through spherical infall, i.e., the 
monolithic collapse of a density fiuctuation of isotropically 
distributed, cold, coUisionless matter in an otherwise ho- 
mogeneous expanding universe. This allowed them to de- 
rive the density profile resulting from self-similar initial 
conditions in an Einstein-de Sitter universe up to the on- 
set of shell crossing. The effects of shell crossing and of 
adopting more and more realistic initial conditions were 
addressed and the whole treatment was refined in a series 
of subsequent papers (Gott 1975; Gunn 1977; Fillmore & 
Goldreich 1984; Bertschinger 1985; Hoffman & Shaham 
1985; Ryden & Gunn 1987; Ryden 1988; Zaroubi & Hoff- 
man 1993; Henriksen & Widrow 1999; Lokas 2000; Lokas & 
Hoffman 2000; Del Popolo et al. 2000; Engineer, Kanekar, 
& Padmanabhan 2000; Nusser 2001). 

According to these studies, after the collapse of some ini- 
tial seed, a stationary regime is reached, the so-called sec- 
ondary infall, in which the system grows inside-out. This is 
the consequence of the fact that the orbital period of par- 
ticles in any shell is smaller than the characteristic time of 
secular variation of their apapsis (or turn-around radius) 
owing to shell crossing of the infalling layers. Accordingly, 
1) particles spend most of the time near the apapsis of 
their orbits, which makes their time-averaged radius close 



to that value, and 2) shell crossing gently alters the initial 
apapsis (corresponding to maximum expansion) of layers, 
making them contract by a monotonously varying factor of 
order two and, hence, with no apapsis crossing. This is an 
important result since, provided the infall rate of matter 
is known, the inside-out growth condition completely de- 
termines the halo density profile. Much progress has been 
achieved in the last twenty years in the modeling of halo 
mass growth. In particular, the extended Press-Schechter 
(PS) model (Press & Schechter 1974; Bower 1991; Bond 
et al. 1991; Lacey & Cole 1993) makes accurate predic- 
tions on the rate at which halos increase their mass in 
hierarchical cosmologies (Lacey & Cole 1994). Unfortu- 
nately, in such cosmologies, halos develop from small to 
large scales by successive aggregations rather than through 
smooth spherical infall. 

There are different ways in the literature to refer to the 
mass assembly process. Some authors call each individ- 
ual aggregation event "merger" (e.g.. Bower 1991; Bond 
et al. 1991; Lacey & Cole 1993), while others use the word 
"accretion" to refer to their cumulative effect (e.g., Wech- 
sler et al. 2002, hereafter W02; Zhao et al. 2003). In N- 
body or Monte Carlo simulations, it is also usual to use 
the words "merger" and "accretion" depending on whether 
the captured clumps are resolved or not, respectively (e.g., 
Somerville & Kolatt 1999). On the other hand, if one is 
only interested in the mass growth of halos, the distinc- 
tion between "major" and "minor" mergers or "fast" and 
"slow" accretion is not important as it represents solely an 
arbitrary classification of captures according to the mass 
increase they produce. However, in dealing with the struc- 
tural and kinematical effects of aggregation events on a 
given halo, such a distinction is crucial. 

When a halo captures a much less massive clump its 
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structure is essentially unperturbed while, in the case 
that both objects have comparable masses, the system is 
brought far from equilibrium and violently relaxes again 
(Lynden-Bell 1967), thus loosing the memory of the pre- 
vious history. The new equilibrium state is characterized 
by a normal distribution of the velocities of the differ- 
ent constituents similar to that yielded by two-body re- 
laxation, although independent of the particle mass. Con- 
sequently, were halos unperturbed after any such dramatic 
event, they would end up as spherical systems with uni- 
form, isotropic velocity dispersion. For this reason, halos 
are often modeled as isotropic, monomassic, isothermal 
spheres (e.g., King 1972; Shapiro, Iliev, & Raga 1999). 
However, halos are not isolated systems and they continue 
collecting matter during violent relaxation, which makes 
the final state difficult to predict (e.g., Saslaw 1987). In 
addition, the frontier between minor and major captures 
is rather blurred and intermediate captures have their own 
specific dynamical effects. The associated clumps undergo 
substantial dynamical friction and spiral to the halo cen- 
ter rather than move along stable orbits. Then, the den- 
sity distribution is neither unaltered nor completely rear- 
ranged, but becomes a little cuspicr at the halo center. 

The situation is therefore quite complicated and diffi- 
cult to implement in an analytical model. For this rea- 
son some authors have turned to numerical experiments. 
Navarro, Frenk, & White (1996; 1997, hereafter NFW) 
found, through A^-body simulations of several hierarchical 
cosmologies, that the spherically averaged density profile 
of current halos of a wide range of masses is always well 
fitted by the simple analytical expression 
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hereafter referred to as the NFW profile. In equation (1), 
r is the radial distance to the halo center, and Tg and pc 
are the halo scale radius and characteristic density, re- 
spectively. The latter two parameters are related to each 
other and to the mass M of the halo through the condition 
that the average density within the virial radius R of the 
system is equal to / times pu, with / a constant factor 
independent of mass and some characteristic density of 
the universe. NFW adopted pu equal to the critical density 
Pcrit, and / = 200, although p^ equal to the mean cosmic 
density p and other values of / are also often adopted. 

On the other hand. Cole & Lacey (1997) showed that 
the spherically averaged, locally isotropized, velocity dis- 
persion profile S(r) is well fitted by the solution of the 
Jeans equation for hydrostatic equilibrium and negligible 
rotation 
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where G is the gravitational constant, p{r) the density 
profile given in equation (1) and M(r) the corresponding 
mass profile, for the boundary condition of null pressure 
at infinity. 

These results have been extensively confirmed (e.g., Tor- 
men, Bouchet, & White 1997; Huss, Jain, & Steinmetz 
1999b; Bullock ct al. 2001). There is only some controversy 
on the exact slope of the density profile at very small radii 
(Moore et al. 1998; Jing & Suto 2000; Fukushige & Makino 
2001; Klypin et al. 2001; Power et al. 2003), the most af- 



fected by the limited resolution of simulations. Nonethe- 
less, in spite of the substantial progress in our knowledge 
on the inner structure of halos afforded by iV-body simu- 
lations, the physical origin of the halo density and velocity 
dispersion profiles remains poorly understood. 

Another well known result inferred from iV-body sim- 
ulations is that, in all the cosmologies investigated, the 
halo concentration c = R/rg decreases with increasing halo 
mass. NFW proposed that this is a consequence of the 
fact that, in hierarchical cosmologies, more massive halos 
form later when the mean density of the universe is lower. 
They demonstrated that the empirical mass-density rela- 
tion is automatically recovered when the proportionality 
Pc oc p{t[) is assumed. They were forced, however, to 
define the formation time t{ as the time some progenitor 
collected one hundredth of the final halo mass. Salvador- 
Sole, Solanes, & Manrique (1998, hereafter SSM) showed, 
nonetheless, the robustness of that result using a better- 
motivated definition of the halo formation time: the last 
time the system was rearranged in a major merger. They 
also showed that the relation pc oc p{ti) (strictly with 
Pcrit instead of p, though this makes no difference for the 
typical halo formation times in any popular cosmology) 
is equivalent to a universal dimensionless density profile 
for newborn halos, thus bringing some physical motiva- 
tion to the assumed proportionality. Unfortunately, the 
^-dependence of the concentration predicted by such an 
assumption is not corroborated by A''-body simulations 
(Bullock et al. 2001, hereafter BOl). 

A more promising result found by Avila- Reese, Firmani, 
& Hernandez (1998) was that the density profile of halos 
that have not been subjected to any major merger is very 
similar to the NFW profile. As discussed below, the con- 
tinuous aggregation of small clumps proceeds essentially 
as secondary infall. It is therefore well understood why, 
using the spherical infall model from quite a realistic den- 
sity profile of the initial perturbation such as the typical 
density run around peaks, Del Popolo et al. (2000) also 
found a density profile similar to the NFW one. On the 
other hand, since halos grow inside-out during secondary 
infall, these results are also consistent with the finding by 
Nusser & Sheth (1999) and KuU (2000) that the density 
profile of halos growing inside-out (or through stable clus- 
tering as the former authors call this process) and having 
not experienced a major merger is similar to the NFW 
profile. 

However, the density profiles obtained from spherical 
infall do not satisfy the correct mass-concentration rela- 
tion: the predicted concentration is substantially higher 
than that of simulated halos. Therefore, the role of spher- 
ical infall in the density profile of halos is still debatable. 
In this one respect, is also important mentioning that the 
inside-out growth condition is inconsistent with the dy- 
namical effects of intermediate mass captures as mid-size 
clumps tend to migrate to the halo center. Syer & White 
(1998) have examined the possibility that the NFW profile 
results from the repeated action of this kind of captures, 
concluding that it might explain the central cusp (see also 
Subramanian, Cen, & Ostriker 2000). But, as shown by 
Nusser & Sheth (1999), when the specific effects of in- 
termediate captures are taken into account, the resulting 
density profile deviates more from the NFW profile than 
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the one obtained under strict accretion the role, in this 
scheme, of the tidal stripping of captured clumps has been 
also examined by Dekel, Devor, & Hetzroni (2003); Dekel 
et al. (2003). Thus, before spherical infall can be confirmed 
as the main responsible for the density profile of halos, it 
is necessary to understand why the alterations induced by 
intermediate captures do not show up more frequently. 

More importantly, all these works do not take into ac- 
count the dramatic effects of major mergers (or fast accre- 
tion). By definition, this process causes the rearrangement 
of the system destroying in this way the density profile de- 
veloped until that moment by spherical infall. Of course, 
this would go unnoticed if the density profile emerging 
from violent relaxation turns out to be very similar to 
the one a halo of the same mass would have developed 
through spherical infall. Surprisingly enough, numerical 
simulations indicate that this is just what happens. In- 
deed, high-resolution A^-body simulations show that "the 
occurrence of a recent [major] merger is not an impor- 
tant factor affecting the [shape and the] concentration" 
of simulated halos (W02, p. 568). Likewise, the phase 
space density, p{r)/T,'^{r), of halos in simulations of hi- 
erarchical cosmologies is found to follow, independently 
of their assembly history, a power-law for more than two 
decades in radius with logarithmic slope close to that pre- 
dicted by Bertschinger (1985) using the spherical infall 
model (Taylor & Navarro 2001). Furthermore, the den- 
sity profile of halos obtained in numerical simulations of 
monolithic collapse are very similar to those obtained in 
hierarchical cosmologies, which indicates that "the merger 
history does not play a role in determining the halo struc- 
ture" (Moore et al. 1999, p. 1147; see also Huss, Jain, & 
Steinmetz 1999a). In other words, the density profile of 
simulated halos is always the same regardless of whether 
(and when) their structure has been rearranged by major 
mergers. This necessarily implies that the density profile 
emerging from violent relaxation must be very similar to 
that developed through spherical infall. But what is the 
reason for such a coincidence? 

The present paper attempts to answer all these ques- 
tions. We begin in § 2 by describing the basic equations 
of a variant of the extended PS model that are used in the 
computations of subsequent sections. We then examine, in 
§ 3, whether pure accretion, which yields density profiles 
of the NFW type, can recover the mass-concentration re- 
lation of simulated halos. In § 4, we analyze the effects of 
mergers and intermediate captures. The results are sum- 
marized and discussed in § 5. 

2. MASS GROWTH THROUGH ACCRETION AND MERGERS 

The Modified Press-Schechter (MPS) model (SSM; 
Raig, Gonzalez- Casado, & Salvador-Sole 1998, and 2001, 
hereafter RGS) is a variant of the extended PS model 
intended to describe not only the mass growth of dark- 
matter halos in hierarchical cosmologies, but their inner 
structure as well. To this end it distinguishes between ma- 
jor and minor aggregation events (i.e., major and minor 
mergers or fast and slow accretion) according to whether 
they cause or not the complete rearrangement of the sys- 
tem, based on the usual comparison of the resulting frac- 
tional mass increase with respect to the reference halo with 
some pre-established threshold. Am, separating the two 



regimes. Note that a given aggregation event may be seen 
as major or minor depending on the partner halo consid- 
ered. For consistency, we say that halos are destroyed 
in major aggregation events while they survive in minor 
ones. We will follow here the most usual notation adopted 
in observational studies (and in SSM) and refer to ma- 
jor aggregation events as "mergers" and to the continuous 
aggregation of small clumps as "accretion" (one should 
avoid confusion with the different meaning of these words 
in other works; § 1). 

According to the definition of Am, the accretion rate, 
that is, the rate at which halos with M at t increase their 
mass between two mergers, is 

r''{M,t)= / (M'-M)rLc(M,M',i)dM', (3) 

J M 

with ri,Q{M, M', t) AM' the instantaneous transition rate 
at t from halos with M to halos between M' and Af'+dAf', 
provided by the extended PS model (Lacey & Cole 1993) 
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In equation (4), 5c{t) is the linear extrapolation to the 
present time of the critical overdensity of primordial 
fluctuations collapsing at t, and a{M) = a{M,to) is the 
r.m.s. fluctuation of the density field at to smoothed over 
spheres of mass M. Both 5c{t) and a{M) depend on the 
cosmology. The M{t) track followed, during accretion, by 
halos with Mj at ti is therefore the solution of the differ- 
ential equation 

^=r^[M(t),t] (5) 

for the initial condition M{ti) = Mj. Strictly speaking, 

this is the average track followed by those halos. Real 
accretion tracks actually diffuse from it owing to the ran- 
dom effects of individual captures. The scatter remains 
nonetheless quite limited along the typical lifetime of ha- 
los (see RGS). 

The destruction of a halo does not necessarily imply the 
formation of a new one: the largest partner participating 
in the merger can perceive it as minor and survive. Only 
those mergers in which all initial halos are destroyed or, 
equivalently, none of them survives do mark the formation 
of a new halo. (The definitions of halo destruction and 
survival imply, indeed, that there is at most one surviving 
halo in any aggregation event.) Thus, the formation of 
a halo corresponds to the last time the system was rear- 
ranged in a merger. 

To derive the formation rate we need to introduce the 
capture rate, r^{M' , M, t) dM, giving the rate at which a 
halo with mass M' at t results from a major merger of a 
halo with mass M to M -I- dM 

r^iM', M, t) dM At = tlc (M, M', t) e[M' -M{1 + Am)] 

x^^dMdt, (6) 
N{M',t) ^ ' 

with rLc(A^,M',t) given by equation (4), N{M,t) the 
PS mass function, and 9{x) the Heaviside function. As 
shown by RGS, mergers leading to the formation of new 
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halos arc essentially binary. This is consistent with the 
fact that the capture rate is closely symmetrical around 
M = M' /2 (see Fig. 1 in RGS), at least in the range 
[M'Am/(l + Am),M7(l + Am)] Corresponding to captures 
leading to the formation of new halos. Note that a frac- 
tional mass relative to the initial object equal to Am rep- 
resents a fractional mass relative to the final object equal 
to Ai„/(1 + Am). Since each symmetric pair of such cap- 
tures produces the formation of one new halo with mass 
M', the formation rate of halos of this mass is 

1 /-M'/Cl+A™) ^^ 
^ Jm'A„,/(1+A„) -'^(M ,t) 

(7) 

The probability distribution function (PDF) of forma- 
tion times can then be derived by taking into account that 
the mass evolution of a halo since its formation is given by 
the accretion track M{t) solution of equation (5). Thus, 
the cumulative number density of halos at ti with masses 
in the arbitrarily small range Mj to Mi-\-5Mi that pre-exist 
at a time t < ti ov, equivalently, the spatial number den- 
sity of halos that evolve by accretion from t to U ending 
up with a mass between Mi and Mj -|- 5Mi is 



7Vp,e(t) = N{Mu ti) SMi exp i - / r^[M{t'), t')] dt' 
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(8) 

Consequently, the PDF of formation times for halos with 
masses between Mi and Mi -|- SMi at ti is given by 



Mt\Mi) 
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N{Mi,ti)SMi dt 
= r^[M{t),t\ exp 



r^[M{t'),t']dt'Y 



(9) 



The median for this distribution defines the typical halo 
formation time. 

The validity of all preceding analytical expressions was 
checked against numerical simulations in RGS. The fact 
that the simulations performed by Lacey & Cole (1994) 
used in that comparison have low resolution by today stan- 
dards should not affect the conclusions drawn in RGS as 
they refer to the way halos grow and not to their inner 
structure. However, it is worth re-examining, now with 
the help of high-resolution simulations, the correct behav- 
ior of the mass accretion rate (eq. [3]) determining the 
density profile of halos in the present model (see next sec- 
tion). 

In Figure 1, we compare the accretion histories predicted 
from the MPS model with those obtained by W02 from the 
high-resolution simulations of BOl. The accretion histories 
of W02 have been converted to the definition of halo mass 
used in the present work (see § 3 for the values of / and Pu 
adopted). Note also that since the latter accretion histo- 
ries include both major and minor aggregations they must 
be compared with our theoretical predictions drawn from 
equation (3) for Am equal to unity. As can be see, there 
is good agreement between theory and simulations. Only 
for halo masses above the critical mass for collapse, M*, 
there is an increasing departure. This seems to reflect the 
difhculty of the PS formalism to correctly model the mass 
growth of halos towards large masses (see, e.g., Monaco 
1997), although no such departure is detected in Fig. 3 of 
RGS. 



3. THE DENSITY PROFILE SET BY ACCRETION 

The rearrangement, through violent relaxation, of a halo 
at its formation produces a more or less spherical system 
with a new density profile independent of the past history 
of the halo. On the other hand, accretion consists of fre- 
quent multiple captures of very small halos. Thus, the 
graininess of accreted matter can be neglected in a first 
approximation. The resulting configuration (i.e., a central 
relaxed spherical object surrounded by a rather smooth 
distribution of matter falling into it) approximately satis- 
fies the conditions met in secondary infall. Consequently, 
halos grow inside-out during accretion, making it possible 
to determine the growth of their density profile from the 
accretion rate derived in the preceding section. 

3.1. Predictions at z = 

Consider a halo with current mass M having formed 
at tf. During the subsequent accretion phase, the halo 
follows the M{t) track solution of the differential equa- 
tion (5). According to the inside-out growth condition, 
the mass accreted at any moment t is deposited at the 
instantaneous radius R{t) of the system without altering 
the inner density profile. Consequently, we have 



rR(t) 

M{t) - Mi ^ 47rr^p(r) dr 

JRf 



(10) 



with Mf and i?f the mass and radius at formation and 
p{r) the density profile developed since that moment. By 
difi^erentiating equation (10) and taking into account the 
definition of the virial radius 



m = 



3M{t) 



/47rpu(t) 



1/3 



and equation (5), we are led to the expression 



M{t) dlnpu 



p{t) = fPnit) < 1 



r'^[M{t),t] dt 



(11) 



(12) 



giving the density at r = R{t). Equations (11) and (12) 
therefore define, in the parametric form, the density pro- 
file developing by accretion from the initial seed at t{, at 
any radius r > Rf. Notice that, owing to the abovemen- 
tioned dispersion of real M{t) tracks around the solution 
of equation (5), the previous theoretical profile represents, 
like the NFW profile, the average density profile of halos 
with M at t, with individual profiles fiuctuating around it. 

The shape of this theoretical profile depends on the val- 
ues of / and pu entering the definition of the virial radius 
(eq. [11]). The need to choose appropriate values for these 
two quantities in order to obtain good predictions from 
the PS model is not new. For instance, the conditional as 
well as the unconditioned mass functions of halos show- 
ing a given abundance depend on the radius adopted to 
define the halo mass. In Figure 2, we plot the theoreti- 
cal density profiles predicted from p^ = pcrit and p in a 
flat, fim = 0.25 CDM cosmology assuming that accretion 
operates from a very small seed. As can be seen, p^ = p 
leads to a monotonous decreasing density profile, while 
Pu = Pcrit leads to a density profile which tends to level off 
at those radii corresponding to cosmic times when A be- 
comes dynamically dominant. Since such radial behavior 
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is not observed in real or simulated halos, we conclude that 
the best value of to use is p. The resulting profile can 
then be further truncated, of course, at any smaller radius 
if necessary. This will be done when comparing our theo- 
retical profile with the NFW profile that uses = pcrii- 
Concerning /, the theoretical density profile is little sensi- 
tive to it, any usual value yielding similar results. Accord- 
ingly, unless we state otherwise, we will use / = 200 as in 
NFW. 

The theoretical density profile also depends (through cq. 
[3]) on the value of Am, the effective threshold for merg- 
ers. As explained in RGS, the comparison between theory 
and simulations in terms of halo masses alone puts no con- 
straint on this parameter. In contrast, the inner structure 
of halos can be used to fix the value of this parameter. In 
Figure 3, we show the effects of varying Am within the 
range [0.3,0.7] around the best value of 0.5 (see below), 
which we use hereafter. 

Once the values of Am, Pu, and / have been fixed, the 
theoretical density profile of halos with M at t having 
evolved by accretion from the initial seed formed at the 
last major merger can be inferred from equations (11) and 
(12). In Figure 4, we show the profiles obtained assuming 
pure accretion^ meaning that the last major merger lead- 
ing to the formation of the halo took place long ago, and 
we compare them to the corresponding NFW profiles. We 
plot only three cosmologies, although similar results are 
obtained in any of the cosmologies studied by NFW. Taken 
as a whole, there is a good agreement between the theo- 
retical and empirical profiles down to r » 10~^i? for halos 
spanning at least four decades in mass. However, a sys- 
tematic tendency is observed for the accordance to worsen 
beyond M*, becoming substantial at about 10 times M*. 
This departure is likely caused by the abovementioned de- 
viation for large halo masses between the accretion his- 
tories predicted by the MPS model and those obtained 
in high-resolution simulations. But the most massive ha- 
los are scarce and seldom virialized (typically they have 
formed recently), so the disagreement between our model 
and the data is not important in practice. 

Both the theoretical and empirical profiles shown in Fig- 
ure 4 for different halo masses are completely fixed in shape 
as well as in scale. This means that our approach yields 
not only a density profile similar to the NFW one (at least 
for moderate and small halo masses), but also good val- 
ues of c. In this respect, it is worth noting that the NFW 
expression is a fit, i.e., only a good approximation to the 
real density profile of simulated halos, and that the em- 
pirical mass-concentration relation is inferred from a small 
number of halos of different masses in each cosmology (the 
few points in the mass-concentration diagram in Fig. 5). 
This affects, of course, the comparison between the theo- 
retical and empirical halo profiles performed in Figure 4. 
Hence, to better assess the quality of the predictions of our 
formalism it is preferable to adjust, through minimiza- 
tion, the theoretical profiles corresponding to different halo 
masses by the NFW expression, as if they were the profiles 
of individual simulated halos, and compare the resulting 
theoretical mass-concentration relation with iV-body data. 

As shown in Figure 5, the theoretical mass- 
concentration relations above ~ 10 M* deviate strongly 
from the empirical data, particularly for power-law spec- 



tra. As indicated by the values, this departure is 
due to the progressive disagreement between the theoret- 
ical profile and the NFW expression with increasing halo 
mass. In contrast, below ~ 10 M*, the empirical mass- 
concentration relationships for the cosmologies studied by 
NFW are well bracketed for Am between 0.3 and 0.7, with 
the lower and upper bounds tending to yield too large and 
too small concentrations, respectively. Hence, were the 
effective frontier between accretion and mergers lowered 
to exclude intermediate captures, the predicted concentra- 
tion would be larger than that found in simulated halos, in 
agreement with the results of Nusser & Sheth (1999), Kull 
(2000) and Del Popolo et al. (2000). In § 4 we wiU provide, 
however, justification for the classification of intermediate 
captures as accretion. 

3.2. Other redshifts 

Since our model can be applied to any redshift, we can 
also derive the redshift dependence of c for halos of a fixed 
mass and compare it to that found by BOl in iV-body sim- 
ulations of the ACDM cosmology. These authors adopted 
a halo radius R defined according to equation (11) with 
Pu — p SiS here, but / varying from 337 to 178 as z shifts 
from to 5 according to the top-hat collapse model in 
the ACDM cosmology. Therefore, in order to carry out 
this comparison one can derive the theoretical profile us- 
ing / = 200 as above, and then truncate or extend it to the 
same z-dependent radius as in BOl or, alternatively, one 
can take advantage of the fact that the theoretical profile 
is little sensitive to / and derive it directly using the same 
f{z) of BOl. In the former case, the halo mass can only 
be known a posteriori after the appropriate truncation or 
extension of the halo density profile, while, in the latter 
case, it coincides with the mass used by BOl. For this 
reason, we have chosen to follow the latter approach. 

In Figure 6 we show the predicted c{z) relation and com- 
pare it to that found by BOl for a halo mass equal to 
5.55 X 10^^ M0, the central logarithmic value of the mass 
range 2.14 x 10^^ - 1.43 x lO" Mq studied by these au- 
thors. This mass is small enough, along the whole range 
of redshifts involved in the comparison, for the theoretical 
profile to always be well adjusted by the NFW expression 
while, at the same time, it is sufficiently large to guaran- 
tee that simulated halos have a large number of particles 
and, hence, that their concentration is well determined. 
The agreement between the predicted and empirical c(z) 
relations is good. We want to stress that this is the first 
correct prediction of c{z) made by a physical model. Pre- 
vious predictions based on the proportionality pc cx Pu(tf) 
proposed by NFW (corresponding to a universal dimen- 
sionless density profile of halos at formation; SSM) have 
been proved, on the contrary, not correct regardless of the 
exact definition of tf used (see Fig. 6), while other analyt- 
ical expressions which correctly fit the empirical function 
c(z) are but toy models (BOl; Eke, Navarro, & Steinmetz 
2001). 

What is wrong then with the assumption pc oc /Ou(if)? 
According to the present model, all halos lying, at differ- 
ent epochs, along a given accretion track M{t) have the 
same density profile, though truncated at their respective 
virial radii. The fact that the theoretical density profile of 
halos at 2: = is reasonably well fitted by the NFW ex- 
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pression over the whole radial range implies that all those 
halos have essentially the same values of the scale param- 
eters Pc and Ts (though not of c, which varies with z owing 
to the time dependence of R). In other words, pc and Vg 
do not change along accretion tracks. Consequently, they 
remain constant for halos evolving by accretion (in agree- 
ment with the results of A''-body simulations; Zhao et al. 
2003). In particular, the value of pc does not depend on 
the formation time of halos, but on their (last) accretion 
track. 

4. THE EFFECTS OF MERGERS AND INTERMEDIATE 

CAPTURES 

The preceding theoretical profile has been derived as- 
suming pure accretion and the inside-out growth of ha- 
los during such a process. This presumes: 1) that halos 
have not undergone any merger for a very long period of 
time and 2) that the capture of clumps with intermedi- 
ate masses has a negligible contribution in the accretion 
process. But both conditions are far from realistic in hi- 
erarchical cosmologies. For instance, A^-body simulations 
(Tormen 1997) show that the average numbers of aggrega- 
tions producing, in the standard CDM cosmology, a frac- 
tional mass increase in halos of cluster scale above 0.1, 
0.2, and 0.5 along the age of the miiverse are, respectively, 
13.7 ± 1.4, 7.4 ± 1.0, and 2.3 ± 0.5. These values reduce 
to 2.7 ± 0.5, 1.5 ± 0.4, or 0.1 ± 0.1 if we only consider the 
time elapsed since the formation of the halos, defined as 
in the original extended PS model, i.e., the time that half 
the mass of the halo was first collected in some progenitor 
(Lacey & Cole 1993), while om- definition of halo formation 
time (see § 2) leads to intermediate values. 

4.1. Mergers 

Though rare, mergers are actually frequent enough to 
invalidate the assumption of pure accretion. This is illus- 
trated in Figure 4. Arrows mark the radii of the seeds 
formed at the typical last merger of halos inside which the 
density profile predicted assuming pure accretion should 
be replaced by that yielded by violent relaxation. How- 
ever, as mentioned in § 1, the latter profile is eventually 
very similar to the one the system would have developed 
through pure spherical infall. We want to emphasize that 
this is not a consequence of our modeling but an empir- 
ical fact drawn from high-resolution cosmological simula- 
tions. As shown next, such a "coincidence" emanates from 
the fact that violent relaxation takes some finite time to 
proceed during which accretion keeps going on. The fi- 
nal density profile then adopts the unique inner structure 
compatible with the boundary conditions imposed by the 
accreting layers. 

We begin by showing first that the velocity dispersion 
profile E(r) of halos in the outer accreted region is uniquely 
determined. It could be derived, of course, from the Jeans 
equation (2) provided we knew the value it takes at some 
given radius. But this is information that we do not have 
a priori (we are avoiding any unjustified assumption at 
infinity). What we only know is that the velocity disper- 
sion profile of seeds emerging from violent relaxation in 
real non-isolated systems deviates from the uniform one 
expected for ideal isolated systems owing to the boundary 
conditions imposed by the infalling matter. This suggests 



the following iterative way to determine the value of S at 
r = the frontier between the inner seed and the outer 
accreted region. 

The zero-order solution of at r = R[ is given by 
the squared value of the uniform velocity dispersion that 
violent relaxation tends to establish, given by equation 
(2) with null logarithmic derivative of and logarithmic 
derivative of p corresponding to the known outer density 
profile. But violent relaxation takes a finite time to pro- 
ceed, during which the instantaneous radius of the; halo 
shifts outwards. Consequently, the same condition will 
hold at the new edge of the system after some arbitrar- 
ily small time. We can therefore derive in the same way 
the zero-order value of at this new edge of the system. 
Those two values of can then be used to determine 
the zero-order logarithmic derivative of at r = i?f and, 
by substituting this derivative into the equation (2), to 
infer the first-order value of at that point. This iter- 
ative procedure rapidly converges to the wanted value of 
S2(i?f). 

Once we know the value of T,{Ri), we can infer the ve- 
locity dispersion profile in the outer accreted region where 
the density profile is known (eqs. [11] and [12]) from the 
Jeans equation (eq. [2]). But this velocity dispersion pro- 
file must coincide with the one that can be obtained by 
directly applying the previous iterative procedure at any 
point of the outer region as both solutions satisfy, by con- 
struction, the Jeans equation for the same density profile 
and have identical values at R{. To sum up, thanks to the 
uniform velocity dispersion profile that violent relaxation 
tends to establish in the inner seed, the velocity disper- 
sion profile emerging in the outer accreted region is fully 
determined, turning out to be also as close to uniform as 
allowed by the density run there. In Figure 7, we compare, 
for the same halo masses and hierarchical cosmologies as 
in Figure 4, the velocity dispersion profile of simulated ha- 
los (Cole & Lacey 1997) with the theoretical one inferred 
following the previous iterative procedure over the whole 
radial range as would correspond to the case of pure ac- 
cretion. As can be seen, the agreement between the theo- 
retical and empirical velocity dispersion profiles is as good 
as in the case of the density profiles. 

We can now proceed with our reasoning and focus on the 
density profile of the inner seed. Since violent relaxation 
lasts for some time during which matter is continuously 
falling into the halo, the final mass distribution adapts 
to the boundary conditions imposed by this matter in- 
fall. Thus, it must be in steady state, compatible with the 
shell crossing of infalling layers, and with a velocity dis- 
persion as close to uniform as possible. But these are the 
same conditions satisfied by the halo during the accretion 
process. It is therefore well understood why the density 
profile emerging from a major merger coincides with the 
profile that would have developed by means of pure ac- 
cretion. Furthermore, since the concept of orbital decay 
is meaningless during violent relaxation, there is no differ- 
ence between intermediate and small captures during that 
process apart from the distinct mass increase they pro- 
duce. Therefore, intermediate captures must be included 
in the formal accretion rate leading to the same density 
profile as violent relaxation. 

It might be argued that the previous reasoning presumes 
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that violent relaxation proceeds to completion, while the 
relaxation time, equal to, say, three crossing times of the 
system at the time of the major merger, 

irei(if) « 3 « 0.66 [GpiU)]-'/' , (13) 

with Vc = {GM/Ry/"^ the circular velocity of the halo, is 
long enough for that condition not always to be satisfied. 
However, the NFW expression describes the average profile 
of halos having a relaxed appearance. That is, those halos 
observed before violent relaxation has gone to completion 
are not taken into account. 

4.2. Intermediate Captures 

Provided that halos do not undergo intermediate cap- 
tures after the completion of violent relax;ation at their last 
major merger, their density profile will develop inside-out 
as explained in § 3. Strictly, there is no need, in the present 
case, for accretion to include intermediate captures, al- 
though given that the outermost profile is quite insensitive 
to the value of Am used (see Fig. 3) this makes almost no 
difference in practice. If, on the contrary, halos undergo 
intermediate captures in such a late phase, as one would 
naively expect from the larger frequency of intermediate 
captures as compared to major ones, the density profile 
will become cuspicr and deviate from the NFW profile. 

What is therefore crucial to understand the shape of 
the NFW profile of relaxed halos is to see that such in- 
termediate captures are quite improbable. This is not in 
contradistinction with the relatively high average numbers 
quoted by Tormen (1997). Here we must only consider 
those intermediate captures restricted to occur after the 
completion of the violent relaxation of the halo accompa- 
nying its last merger and early enough for the captured 
clump to have completed its orbital decay by the time the 
halo is observed. This latter bound is necessary, indeed, 
to guarantee that the halo does not show any substruc- 
ture and can then be identified as a relaxed system. The 
low frequency of intermediate captures subject to these 
two constraints seems to be supported by the results of 
A''-body simulations (Ascasibar et al. 2002). 

Let us therefore calculate the probability V{M) that the 
last intermediate merger of a halo with M at t occurred 
after the completicm of violent relaxation and the time at 
which the eventual intermediate merger should take place 
for the merged clump to reach the halo center by the time 
the halo is observed. This probability is given by 

/■to /-Am/CH-Am) 

V{M) = / ^{t{\M)dtf / dA 

/•to— ftrico(^) 

X / Piast(A,t|M)dt, (14) 

with the integrand giving the joint probability that the 
halo was formed from t[ to t[ + dt[ and underwent the last 
intermediate merger with a clump of instantaneous frac- 
tional mass between A and A -|- dA in the infinitesimal 
interval of time around t. This joint probability is simply 
equal to the product of the probabilities <I>(tf |M) dt{, given 
in equation (9), and Piast(A,t|M) dAdt, calculated in the 
Appendix. 

In equation (14), tf + trei(if) is the time of completion 
of violent relaxation after the last major merger at tf, 



with the relaxation time irei given by equation (13), and 
to ~ ifrico(A) is the time at which the intermediate merger 
should occur for the clump to have reached the halo center 
at to- The characteristic time of orbital decay of a clump 
of mass Ms in a circular orbit of radius r is (e.g., Binney 
& Tremaine 1987) 



tdec{r) 



■ r2 



0.43 ln[M(r)/Ms] G 



(15) 



with Vc{r) = [GM{r)/r]^^^ the circular velocity at r. 
Hence, by integrating t^cc <i-r under the approximation of 
a singular isothermal halo (for which the integral is ana- 
lytical) from to the radius of the initial orbit, R{t), at 
the time t of the merger, we have an estimate of the time 
spent by the clump to spiral to the halo center owing to 
dynamical friction 



tiric{t,A) 



1 



v.{t)R{ty 



0.43 ln[M(t)/Ms] 2GMs 



ln(0.3A) 



(16) 



where M{t) and Vc{t) are, respectively, the total mass and 
circular velocity of the halo at the time of the merger. 
Thus, ffrico(A) is given by the previous expression for a 
time t solution of the implicit equation t + tfric(t, A) =to. 
Note that in deriving equation (16), we have assumed, for 
simplicity, that Ms remains constant and equal to about 
30% of the mass of the merged clump after the initial tidal 
stripping, i.e., Mg = 0.3 M{t) A, with A the fractional 
mass increase of the halo relative to the final object. Ac- 
cording to Dekcl, Dcvor, & Hetzroni (2003), the mass of 
clumps is reduced to 30% of its initial value after their or- 
bits have already achieved a substantial decay. But these 
authors assume clumps with a density profile steeper than 
the NFW one at large radii, which means that their clumps 
are more difficult to tidally strip than real ones. On the 
other hand, we expect the real orbits of merged clumps 
to be typically elliptical rather than circular as we are as- 
suming here, which should diminish their time of orbital 
decay. However, this effect should be balanced by the fact 
that, in the elliptical case, clumps fall deeper in the halo 
at their pericenter and, hence, are more severely truncated 
by tides since the beginning. 

The upper bound for A in equation (14) is not A,„ = 0.5 
but Ani/(1-|-Am) = 0.33 since the fractional mass increase 
A is relative to the mass M{t) of the halo after the merger. 
The lower bound Amin is the value of A yielding a charac- 
teristic time of orbital decay for the initial orbit (eq. [15] 
with r = R[t]), equal to to — [tf + trei(tf)], that is, the 
solution of the implicit equation 



0.27A: 



1 



ln(0.3Ami„) ^Gp{io- [if + trci(tf)]} 



to-[tf+trel{tf)]. 

(17) 

Values of A smaller than this limit correspond to clumps 
that do not suffer any significant orbital decay. 

In Figure 8, we show, for different halo masses and cos- 
mologies, the quantity V giving the probability that inter- 
mediate captures can affect the halo density profile that 
would emerge from the last major merger and the subse- 
quent accretion phase, calculated in § 3. As can be seen, 
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this probability decreases with increasing mass. It be- 
comes neghgible for hales more massive than about M* 
in all cosmologies, but even for the smallest halo masses 
considered by NFW, equal to ~ IQ-^M* for CDM spec- 
tra or ^ for power-law ones, it is quite small. 
Certainly, the exact values shown in this plot should not 
be taken too literally given that they depend on the ex- 
act definitions of irci and ifric adopted. But it is equally 
true that equation (14) gives an upper bound for V since, 
in its derivation, it has been implicitly assumed that the 
mass increase in captures comes from one unique clump, 
while it might also be produced by two or more less mas- 
sive clumps not falling then into the category of interme- 
diate captures. What is more important, more than one 
intermediate merger is required to affect appreciably the 
density profile, while the probability that two or more in- 
termediate captures take place in the right time interval is 
obviously much smaller. 

5. SUMMARY AND DISCUSSION 

We have presented a simple analytical model for the 
density profile developed by halos during accretion based 
on their inside-out growth in that regime. The resulting 
density profile is found to be similar to the NFW profile 
in agreement with previous works. However, contrarily to 
these same works, we also find good agreement between 
theory and simulations concerning the mass-concentration 
relations of halos at z = 0. More importantly, our phys- 
ical model correctly predicts the mass-concentration rela- 
tion at any redshift. An analytical model such as the one 
presented here, able to make reasonable predictions on the 
density profile of halos at any epoch and in any cosmology, 
should be then very useful in the modeling of galaxy for- 
mation and evolution. Caution must be taken, however, 
in using the present model on halos more massive than 
~ 10M,(f), with M^{t) the typical mass for collapse at t. 
For masses this large, the theoretical profile shows a sub- 
stantial departure from the empirical one likely due to a 
deficient modeling of halo growth by the PS formalism at 
the large mass end. 

The success of our model compared to previous ones 
appears to rely on the inclusion of intermediate-mass cap- 
tures in the accretion regime despite the fact that they 
do not satisfy the inside-out growth condition. This is 
justified by the realization that the halo density profile is 
actually not set during accretion (except for the narrow 
outermost radial range, in logarithmic units), but at the 
time of the last major merger, both processes leading to 
halos with essentially the same shape. Given that there is 
no difference between the dynamical effects of minor and 
intermediate captures during the violent relaxation result- 
ing from mergers, it is then understandable why interme- 



diate captures must be included formally in the accretion 
process. 

The fact that major mergers yield essentially the same 
density profile than pure accretion is an important re- 
sult of high- resolution numerical simulations that has been 
pointed out only by a few authors (W02; Moore et al. 1999; 
Huss, Jain, & Steinmetz 1999a). We have shown here that 
this unexpected coincidence is simply due to the unicity of 
the steady structure compatible with the boundary con- 
ditions imposed by the ever falling surrounding material. 
Thus, we have been able to reconcile the well-known re- 
sult that pure accretion (or spherical infall) yields density 
profiles of the NFW form with the fact that, in hierarchi- 
cal cosmologies, halos endure from time to time important 
mergers (or short periods of intense accretion) that bring 
them far from equilibrium and cause the rearrangement of 
their preceding structure. Besides, we have demonstrated 
that there is no time for intermediate captures to alter the 
density profile a la NFW established after the last merger 
without giving the halo a non-relaxed appearance. 

Finally, we want to remark that our aim here was not 
to infer an accurate density profile for dark-matter halos, 
but to investigate its possible origin. For this reason we 
have considered the simplified case of pure dark-matter ha- 
los while, in the real universe, about 10 — 15% of the halo 
mass is in the dissipative baryonic component. This might 
have appreciable effects on the density profile of real sys- 
tems. Likewise, we have assumed spherical symmetry and 
neglected halo rotation, as well as the possible anisotropy 
of the local velocity tensor. Real halos have, on the con- 
trary, some angular momentum and are immersed in large 
filamentary structures making them accrete matter pref- 
erentially along one privileged direction (West 1994) and 
feel the tidal field of such anisotropic structures (Salvador- 
Sole & Solanes 1993). Moreover, even with perfect spheri- 
cal symmetry, some velocity anisotropy will prevail in the 
halo outskirts, as observed, owing to the distinct evolution 
of the radial and tangential velocity dispersions of infalling 
layers. We note, however, that our density profile is inde- 
pendent of the actual degree of anisotropy of the velocity 
tensor. 
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APPENDIX 

PROBABILITY OF THE LAST INTERMEDIATE MERGER 

Let us consider the rate of intermediate captures up to some given fractional mass increase A of halos with final mass 
M' at i 

r'iM',<A,t)= r^c{M,M',t)^j^dM. (Al) 

iM'A„i„ I\[M ,t) 

Following the same reasoning leading to the PDF of formation times, but adapted to the previous range of intermediate 
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captures, we have that the cumulative number density of current hales with masses in the arbitrarily small range M to 
M + 5M that have undergone the last intermediate merger up the fractional mass increase A at some time smaller than 
t is 



A^iast(< A, < 1 1 M) = N{M,tQ) SM exp | - ^ r'[M{t'), < A,t')] dt' 



(A2) 



where M[t) stands for the accretion track calculated using a modified threshold for mergers Am corresponding to Amin to 
take into account that, on the final phase, accretion should not include intermediate captures. Hence, the PDF of times 
that halos with current mass M underwent their last intermediate merger yielding a fractional mass increase from A to 
A + dA is simply 



Plast(A,i|M)dA = 



1 



N{M,to)SM 



dtdA 



iVlast(< A,<i|M) 



dA. 



(A3) 
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Fig. 1. — Mass accretion histories for three halos in the flat ACDM cosmology studied by BOX. The MPS predictions for Am equal 1 (thick 
dashed lines) and 0.5 (thin dashed line) arc compared to the analytical expression proposed by W02 to fit the mass accretion histories drawn 
from the Af-body simulations of BOl (full line). M» = 2.39 x 10^"^ Mq is the typical mass for collapse of density fluctuations at z = for this 
cosmology. 




log (r/R) 

Fig. 2. — Theoretical density profiles for halos at z = with masses equal to lO^'^Mt, lO^^M*, M*, and 10 M« (from bottom to top) in a 
flat, Sim = 0.25, CDM model inferred using pu — Pcrit (dashed lines) and pu — p (solid lines). We are using the Scune cosmological parameters 
and delimitation of halos as in NFW and Am = 0.5. M* is the typical mass for collapse of density fluctuations and p* stands for pcrit^* 
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log (r/i?) 

Fig. 3. — Theoretical density profiles for lialos at 2 = with the same masses as Fig. 2 obtained for different values of the threshold for 
mergers: Am = 0.3 (dot-dashed lines), 0.5 (solid lines), and 0.7 (dashed lines). The three cosmologies represented arc a flat, A = 0, CDM 
model (SCDM), a fiat, Qm = 0.25, CDM model (ACDM), and a flat, A = 0, model (FPL) with power-law spectrum of density fluctuations of 
index n = —1. All model parameters are the same as in the corresponding cosmologies of NFW. 




log (r/R) 

Fig. 4. — Theoretical (solid lines) and NFW (dotted lines) density profiles for halos at z = for the same masses used in the previous 
Figures and the same cosmologies depicted in Fig. 3. Arrows mark the halo radii at the typical halo formation times. 
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Fig. 5. — Mass-concentration relation for halos at z = obtained by fitting the theoretical density profiles derived in the present paper 
by the NFW analytical expression. The solid line corresponds to Am = 0.5, while the dot-dashed and dashed lines correspond to the rather 
extreme values of 0.3 and 0.7, respectively. The dark points are those obtained by NFW from high resolution numerical simulations. We 
plot all the cosmologies studied by these authors, which, apart from those mentioned in the previous Figures, are two Clm = 0.10 power-law 
models (OPL) with n = 0, —1, and three more flat, A = 0, power-law models (FPL) with n = 0, —0.5, —1.5. We also show the predictions 
for an open, Qm = 0.25, CDM cosmology (OCDM), normalized as the ACDM model, not included in the study by NFW. 
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Fig. 6. — Concentration vs. redshift in the ACDM cosmology studied by BOl for a fixed halo mass M equal to 5.55 X lO^^ Mq. The 
function proposed by BOl to fit their empirical data (solid line) and the corresponding intrinsic 68% spread (thin solid lines) are compajred to 
our theoretical prediction (dashed line) and the c{z) dependence suggested by NFW (long-dashed line) and SSM (dotted line). 




Fig. 7. — Same as Fig. 4 but for the velocity dispersion profile. 
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Fig. 8. — Probability that currently relaxed halos with mass M have undergone some merger with a clump of mass in the range required 
to modify the NFW-likc profile arising from their last major merger and the subsequent accretion phase. The different curves correspond to 
the cosmologies plotted in Figs. 3, 4, and 7: SCDM (solid line), ACDM (dashed line), and FPL n = -1 (dotted line). 



